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Abstract 

Simple formulas for optimization of VUV and X-ray SASE FELs are presented. 
The FEL gain length and the optimal beta-function are explicitly expressed in 
terms of the electron beam and undulator parameters. The FEL saturation length 
is estimated taking into account quantum diffusion in the undulator. Examples 
of the FEL optimization are given. Parameters of a SASE FEL, operating at the 
Compton wavelength, are suggested. 
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1 Introduction 



Successful operation of the VUV (vacuum ultraviolet) FEL (free electron laser) at the 
TESLA Test Facihty at DESY [1,2], based on SASE (self-amplified spontaneous emission) 
principle [3], has stimulated a rapidly growing interest in the development of VUV and 
X-ray FELs. A number of projects (see, for instance, [4-8]) are now at different stages of 
design and construction. 

At the first stage of a SASE FEL design one looks for the dependence of the FEL saturation 
length on the wavelength, electron beam parameters, undulator parameters, and beta- 
function. Usually the parameters are optimized for the shortest design wavelength since 
the saturation length is the largest in this case. The saturation length is proportional to 
the gain length (e-folding length) of the fundamental transverse mode (see [9] for more 
details). The gain length can be found by the solution of the FEL eigenvalue equation. 

The eigenvalue equation for a high-gain FEL, including diffraction of radiation, emittance, 
and energy spread, was derived in [10,11]. There exist approximate solutions [12,13] of 
this equation. The exact solution was presented in [14] as well as an approximate solution 
(with a limited validity range). The latter solution was fitted [14] using 3 dimcnsionlcss 
groups of parameters, and 19 fitting coefficients. An approximate solution, that fits the 
exact solution in the entire parameter space with high accuracy (better that 1 %), was 
presented in [15]. A numerical algorithm for finding this approximate solution is very fast 
and robust. It was used to obtain the main results of this paper. 

In this paper we present the explicit, simple and rather accurate dependencies of the FEL 
gain length on the beam and undulator parameters. Our formulas are not universal, but 
they provide a good accuracy (better than 5 % for the gain length) in a typical parameter 
range of VUV and X-ray FELs. We present the formulas without derivation since they 
were not derived analytically. In some sense the parametric dependencies were guessed, 
and then the fitting coefficients were found from the solution of the eigenvalue equation. 
For instance, we used only 2 fitting coefficients for the gain length with the optimized beta- 
function. Our formulas allow one to quickly estimate FEL saturation length, including the 
principal effect of energy diffusion in the undulator due to quantum fluctuations of the 
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undulator radiation. In addition, we present two practical examples of using our design 
formulas: optimization of SASE FEL with negligible energy spread, and the limitation on 
SASE FEL wavelength taking into account quantum diffusion. In particular, we suggest for 
the first time the set of parameters for a SASE FEL operating at the Compton wavelength. 



2 Gain length for the optimized beta-function 



Let us consider an axisymmetric electron beam with a current /, and a Gaussian dis- 
tribution in transverse phase space and in energy [14,15]. The focusing structure in the 
undulator is a superposition of the natural undulator focusing and an external alternating- 
gradient focusing. The eigenvalue equation [14,15] is valid under the following condition 
[15]: 

<^ mm 1, 



27r/3 V ' 27rey 

where Lf is the period of the external focusing structure, /3 is an average beta-function, e 
is the rms emittance of the electron beam, and is the FEL resonant wavelength. The 
resonance condition is written as: 

- ^2 ■ (1) 

Here Aw is the undulator period, 7 is relativistic factor, and K is the rms undulator 
parameter: 

K = 0.934 Aw[cm] Brrn.[T] , (2) 



Brms being the rms undulator field. 

In what follows we assume that the beta-function is optimized so that the FEL gain 
length takes the minimal value for given wavelength, beam and undulator parameters. 
Under this condition the solution of the eigenvalue equation for the field gain length ^ can 

^ There is also a notion of the power gain length which is twice shorter. 
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be approximated as follows: 



L,^L,o(l + 5), (3) 



where 



and 



131 ^ .5^ 

I xV'xT {KAjj^ii + K-^y/^ ■ ^ ' 



The following notations are introduced here: = 17 kA is the Alfven current, e„ = 'je is 
the rms normalized emittance, = a^/mc^ is the rms energy spread (in units of the rest 
energy), Ajj = 1 for a helical undulator and Ajj = Jo{K^ /2{1 + K^)) - Ji{K^ /2{1 + K^)) 
for a planar undulator, Jq and Ji are the Bessel functions of the first kind. 

The formula (3) provides an accuracy better than 5 % in the domain of parameters defined 
as follows 

1<^<5 (6) 



A, 

5 < 2.5 ^ 1 - exp ^ ^^^^ 



2 V A 



(7) 



Note that the condition (6) is usually satisfied in reahstic designs of VUV and X-ray FELs 
when one does optimization for the shortest wavelength (defining the total undulator 
length). The condition (7) is practically not a limitation. To illustrate the accuracy of the 
formula (3) we present a numerical example. The following nominal operating point is 
chosen: = 1 nm, A^ = 3 cm, K = 1, I = 2.5 kA, e„ = 2 /im, = 1 MeV, energy is 2.8 
GeV, undulator is planar. We scan over different parameters and compare the gain length 
calculated with formula (3) and by solving the eigenvalue equation [15]. The results are 
presented in Figs. 1-6. 
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We also present here an approximate expression for the optimal beta- function (an accuracy 
is about 10 % in the above mentioned parameter range): 

/J \ p3/2\l/2 

/3opt^ll.2^) ^Vt^(1 + 8^)-'/' (8) 



XrKAjj 

Note that dependence of the gain length on beta- function is rather weak when (5 > /3opt- 

Finally, let us note that the saturation length cannot be directly found from the eigenvalue 
equation. However, with an accuracy 10-20 % one can accept the following estimate: 



^sat 



10 L, (9) 



3 Influence of quantum diffusion in an undulator on saturation length 

Energy spread growth due to the quantum fluctuations of the spontaneous undulator 
radiation can be an important effect [3,16] in future SASE FELs. The rate of the energy 
diffusion is given by [17]: 

^ = ^Aere7^4^^^(^) , (10) 
where Ac = 3.86 x 10"^^ cm, = 2.82 x 10"^^ cm, = 27r/Aw, and 

F{K) = 1A2K + (1 + 1.50K + 0.95K'^)-^ for helical undulator 

F{K) = 1.7QK + {1 + 1.88K + 0.80K^)-^ for planar undulator (11) 

To estimate the FEL saturation length, we accept the following scheme. First, we neglect 
energy diffusion and find a zeroth order approximation to the saturation length from (9), 
(3)-(5). Then we calculate an induced energy spread in the middle of the undulator from 
(10), add it quadratically to the initial energy spread, and find a new expression for 5. 
Then, using (9), (3)-(5), we find the first approximation to the saturation length. Then 
we do the next iteration, etc. Finally, the saturation length can be estimated as 

L,at ^ 10 Lgo , (12) 
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where 

^ = >^ 1" (tJ ^F7J];vi — ^15^ — (13) 

Note that in the latter formula the powers are somewhat simplified. Comparing Eqs. (9) 
and (12), we can introduce an effective parameter 



(14) 



which should be used instead of S in (7) to check the applicability range and in (8) to 
estimate the optimal beta-function. 

Although formula (12) is rather crude estimate, it can be used for quick orientation in 
the parameter space with a posteriori check using a numerical simulation code. 

4 Examples of SASE FEL optimization 

4-1 Optimized FEL with a negligible energy spread 

Formulas, presented in the previous Sections, can be used for the optimization of undulator 
parameters as soon as a specific type of the undulator is chosen. We demonstrate such a 
possibility with the planar NdFeB undulator of which magnetic field can be described by 
the following formula [8]: 

2" 



5max[T] = 3.694 exp 



-5.068^- + 1-52 

V Aw 



for 0.1 < g/X^ < 1 , (15) 



where g is the undulator gap. The rms value of the parameter K is given by Eq. (2) with 

-^rms -^max/ 

We assume that the energy spread effect on the FEL operation can be neglected {S , Sg ^ 
0). Then, using (3), (2) and (15), we minimize the gain length for a given undulator gap. 
The optimal undulator period is found to be 

(Aw)opt[ciii] ~ 1 + 2^1 [cm] for g > 0.5 cm (16) 
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The optimal value of K is then defined from (15) and (2), the electron beam energy - from 
(1), and the optimal beta-function - from (8). The minimal gain length can be expressed 
(in practical units) as follows: 

(L ) . fml ~ 20 4/n/^m]yV^[cm] 

Using estimate of the saturation length (9), we find the minimal wavelength at which 
SASE PEL can saturate within the given undulator length L^: 

(X) • [Al ~ 3 X 10^ 4''1/^m]^3/4[^^] 

(A.)_[AJ _ 3 X 10 ^3/,^j^^j^3/2j^j (18) 
4-2 SASE PEL at the Compton wavelength 



Another example is the optimization of sub-Angstrom FELs for which the effect of quan- 
tum diffusion in the undulator can play an important role. Wc consider the case when the 
energy spread is dominated by the quantum diffusion, and neglect initial energy spread 
{S 0). Optimizing undulator period and parameter K in (12), we get the following 
estimate for the minimal wavelength ^ : 

(A.)^,[A] - '^"[^T/l (19) 

Note that in some cases the optimal undulator parameters can be impractical. In any 
case, the estimate (19) gives a lower hmit. The following numerical examples show that 
one can be close to this limit with technically feasible undulator parameters. 

Let us consider the electron beam parameters (peak current and emittance) assumed in 
[18]. One of the examples, considered in [18], is a SASE PEL operating at = 0.28 A 
with / = 5kA and e„ = 0.3 /im. Another example is even more ambitious: A,. = 0.12 A 
with 1 = 5 kA and e„ = 0.1 /xm. 

We try to push the wavelength closer to the extreme given by Eq. (19). In our first 
example we assume 1 = 5 kA and e„ = 0.3 /im. With these parameters the wavelength 



One can notice the difference with more crude estimate presented in [16] 
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Xr = 0.1 A can be reached at the electron beam energy 23 GeV in a planar undulator 
with Aw = 2 cm and K = 1 (with the gap g = 0.7 cm according to (15) and (2)). The 
optimal beta-function is about 40 m, and the saturation length is estimated at 160 m. 

The second example is a SASE FEL operating at the Compton wavelength, A,. = Ac = 
0.0234 A (photon energy is 0.5 MeV). We assume the electron beam with / = 5kA and 
en = 0.1 fim, the energy is 40 GeV. We choose a helical undulator with Aw = 2 cm and 
K = 0.7. The optimal beta-function is about 35 m, and the saturation is reached within 
200 m. Our estimates show that quantum effects, other than energy diffusion, give small 
corrections to the classical description and can be neglected. 
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Fig. 1. Gain length versus resonant wavelength for the following set of parameters: = 3 cm, 
K = 1, I = 2.5 kA, e„ = 2 |Ltm, = 1 MeV. Undulator is planar, resonance is maintained by 

tuning electron beam energy, beta-function is optimized for each case. Line is the solution of the 
eigenvalue equation [15], and the circles are calculated using formula (3). 
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Fig. 2. Gain length versus undulator period for the following set of parameters: = 1 nm, 
K = 1, I = 2.5 kA, e„ = 2 /xm, cr^ = 1 MeV. Undulator is planar, resonance is maintained by 
tuning electron beam energy, beta-function is optimized for each case. Line is the solution of the 
eigenvalue equation [15], and the circles are calculated using formula (3). 
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Fig. 3. Gain length versus undulator parameter K for the following set of parameters: = 1 nm, 
Aw = 3 cm, / = 2.5 kA, = 2 ^um, = 1 MeV. Undulator is planar, resonance is maintained 

by tuning electron beam energy, beta-function is optimized for each case. Line is the solution of 
the eigenvalue equation [15], and the circles are calculated using formula (3). 
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Fig. 4. Gain length versus normalized emittance for the following set of parameters: A^ = 1 nm, 
Aw = 3 cm, K = 1, I = 2.5 kA, = 1 MeV. Undulator is planar, beta-function is optimized 
for each case. Line is the solution of the eigenvalue equation [15], and the circles are calculated 
using formula (3). 
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Fig. 5. Gain length versus current for the following set of parameters: A,- = 1 nm, = 3 cm, 
if = 1, e„ = 2 fim, CTg = 1 MeV. Undulator is planar, beta-function is optimized for each case. 
Line is the solution of the eigenvalue equation [15], and the circles are calculated using formula 
(3). 
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Fig. 6. Gain length versus energy spread for the following set of parameters: A^ = 1 nm. Aw = 3 
cm, K = 1, I = 2.5 kA, = 2 /xm. Undulator is planar, beta-function is optimized for each 
case. Line is the solution of the eigenvalue equation [15], and the circles are calculated using 
formula (3). 
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